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Abstract 

We compute the corrections to finite-size scaling for the A^-vector model on the 
square lattice in the large- A^ limit. We find that corrections behave as logL/L^. 
For tree-level improved hamiltonians corrections behave as . In general /-loop 
improvement is expected to reduce this behaviour to log' L). We show that 

the finite-size-scaling and the perturbative limit do not commute in the calcula- 
tion of the corrections to finite-size scaling. We present a detailed study of the 
corrections for the i?P°°-model. 



1 Introduction 



In the study of statistical models it is extremely important to understand finite-size 
corrections. Indeed in experiments and in numerical work it is essential to take into 
account the finite size of the system in order to extract correct infinite- volume predictions 
from the data. Finite-size scaling (FSS) 11 S B H concerns the critical behaviour of 
systems in which one or more directions are finite, even though microscopically large, 
and it is therefore essential in the analysis of experimental data in many situations, 
for instance, for films of finite thickness. Numerically FSS can be used in a variety of 
ways to extract informations on infinite-volume systems. A very interesting method to 
extract critical indices comparing data on lattices of different sizes was introduced by 
Nightingale 0, the so-called phenomenological renormalization group. Recently FSS 
has been used to obtain precise predictions at very large values of the correlation length 
from simulations on small lattices. This extrapolation technique was introduced by 
Liischer, Weisz and Wolff [0 and subsequently applied to many different models ^ ^, 
10| , |TT|, [12], |13[: a careful theoretical analysis (see Sect. 5.1.2 of Ref. [|13|) shows that the 



method is extremely convenient for asymptotically free theories and indeed one was able 
to simulate the 0(3) cr-model [0 up to ^ ^ 10^ and the SU{3) chiral model |I2|, |1|] up 
to ^ ~ 4 ■ 10^ using relatively small lattices {L < 512). In order to use these techniques 
reliably it is extremely important to have some theoretical prediction on the behaviour of 
the corrections to FSS. One can use this information in two different ways. A possibility 
is to take advantage of the theoretical prediction to extrapolate the Monte Carlo data 
to the FSS limit — that still involves a limit L ^ oo — in the spirit of Ref. One 
also needs this information if one determines the FSS curve by comparing data from 
simulations on lattices of different sizes as proposed in Ref. p. For instance checking 
the absence (within error bars) of corrections to FSS for lattices of sizes 64 < L < 256 is 
enough if the corrections vanish as 1/L^ while it can be totally misleading if corrections 
behaving as 1/logL are present (this is the case for instance of the four-state Potts 
model, see Ref. [[T^)- 



A second topic that will be extensively discussed in this paper is the improvement of 
lattice hamiltonians [0, |16|, 0, |18|, [1^ . The idea behind all these attempts is to modify 
the lattice hamiltonian with the addition of irrelevant operators in order to reduce lattice 
artifacts: in this way one hopes to have scaling and FSS at shorter correlation lengths. 
For general statistical models this is a non-trivial program. For asymptotically-free 
theories the idea is much simpler to implement since in this case improvement can be 
discussed using perturbation theory. 

In this paper we will study the problem of corrections to FSS and improvement in the 
context of the large- A^-vector model. This theory provides the simplest example for 
the realization of a nonabelian global symmetry. Its two-dimensional version has been 
extensively studied because it shares with four- dimensional gauge theories the property 



of being asymptotically free in the weak-coupling perturbative expansion |2^, [21], |22 
This picture predicts a nonperturbative generation of a mass gap that controls the 
exponential decay at large distance of the correlation functions. 

Besides perturbation theory, the two-dimensional A^-vector model can be studied 
using different techniques. It can be solved in the N = oo limit |2^ and 1/A^ 
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27 . An exact S'-matrix can be 



corrections can be systematically calculated |2^, g6| 
computed ^ and, using the thermodynamic Bethe ansatz, the exact mass-gap of 
the theory in the limit (3 ^ oo has been obtained 
object of extensive numerical work 



34, 



n|. The model has also been the 
5] mainly devoted to checking the 



correctness of the perturbative predictions [^, 0, [38(1 . 

In the large- limit FSS can be studied analytically 40, 41|. Here we will con- 
centrate on the corrections to FSS in two dimensions and we will compute the deviations 
from FSS for generic lattice interactions. We will show that in general corrections to FSS 
behave as log L/L^. This is in agreement with a general renormalization-group argument 



that shows that corrections to FSS are controlled by the first subleading operator [42 



Tree-level improvement changes the behaviour by a logarithm of L: these actions have 
corrections behaving as Subsequent improvement should reduce the corrections 

to l/{L^\ogL) and so on. Full O(a^) improvement to all orders of perturbation theory 
provides an action with corrections behaving as logL/L^. 

Besides the standard A^- vector model we will also discuss a mixed N-vector/ RP^~^- 
model P^ , There are two reasons why we decided to include this computation: 
first of all, for large values of A^, models with only two-spin interactions show many 
simplifying features: for instance in the /3-function only the leading term does not vanish. 
For this reason one may expect that the behaviour of the corrections for this class of 
models is far simpler than in generic models. Instead the mixed A^-vector/i?P^~^-model 
shows a more complex behaviour and, for instance, the /3-function is non trivial to all 
orders of perturbation theory. We find that in the mixed model the corrections behave 
as (logL/L^)/(L) where f{L) is a non-trivial function such that f{oo) is finite and 
that admits an asymptotic expansion in powers of 1/logL. The presence of powers 
1 / (L^ log" L) is somewhat unexpected from the point of view of perturbation theory. 
We will show that this is related to the non commutativity of the limits L —>■ oo and 
(3 —>■ oo. In other words the perturbative limit /5 — > cxd at L fixed followed by the limit 
L — > cxo gives results that are different from the FSS limit. The commutativity of these 
two limits has been the object of an intense debate. The standard wisdom is that the 
two limits are identical, but this point has been seriously questioned by Patrascioiu and 
Seller (for an answer to these criticisms see Ref. together with many other 

assumptions derived from perturbation theory |^6[. Our calculation shows that, if the 
standard assumption is true, it is a result far from obvious: indeed the limits are different 
for the corrections to FSS. 

on the RP^~^ models that has 



A second motivation lies in recent work |47, 48 



shown the possibility that these models have very large FSS corrections. We wanted to 
understand if there is any sign of this phenomenon in the large- A^ limit. Our explicit 
calculation shows that i?P°° has corrections that are larger than those of the A^-vector 
model. Depending on the observable, for reasonable lattice sizes, we find an increase by 
a factor of 6-15. 



47|. 



This is in qualitative agreement with the scenario of Ref. 
The paper is organized as follows: in Sect. 2 we define the models we consider, and 
compute various observables in the large- A^ limit. In Sect. 3 we discuss the corrections 
to FSS for the A^-vector model and in Sect. 4 we extend our results to the mixed 
N-vectoT/RP'^~^ model. In Sect. 5 we present our conclusions. 

In the Appendices we report some general results on the FSS behaviour of lattice 
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sums. These results are of general interest and may be applied in many other contexts: in 
particular they may be used to study FSS properties of models that have a height (SOS) 
representation (see Refs. p9| , Q and references therein). In Appendix A we define a set 



of basic functions that appear in all our results and we report some of their properties. 
We extend here the results of Ref . |^ . In Appendix B we give an algorithmic procedure 
that allows to compute the expansion in powers of of any sum involving powers of 
the lattice propagator for a Gaussian model with arbitrary interaction in the FSS limit. 
As an example we report the explicit formulae that are needed in our main discussion. 
In Appendix C we report the asymptotic behaviour of some lattice integrals. 



Preliminary results of this work were presented at the Lattice '96 Conference [52 



2 The models 

In this paper we will study the finite-size-scaling properties of the classical A^-vector 
model on a square lattice with local, translation- and parity-invariant ferromagnetic 
interactions. The hamiltonian is given by 

n= NY,J{x-y)(T,-(Ty , (2.1) 

x,y 

where the fields (t^ satisfy cr^ = 1. The partition function is simply 

Z= fYlda^e-^-^. (2.2) 



We will consider general local interactions. If J{p) is the Fourier transform of J{x), 
locality and parity invariance imply that J{p) is a continuous function of p, even under 
p — > —p. We will require invariance under rotations of 7r/2, that is we will assume J{p) 
symmetric under interchange of pi and p2 ■ Redefining f3 we can normalize the couplings 
so that 

J(g) = J(0)-^ + O(g^) (2.3) 

for g ^ 0. We also introduce the function 

wiq) = -2(J(g) - J(0)), (2.4) 

that behaves as q"^ for g — 0. Finally we will require the theory to have the usual 
(formal) continuum limit: we will assume that the equation w{q) = has only one 
solution for —n < qi < tt, namely g = 0. We will need the small-g behaviour of w{q): 
we will assume in this limit the form 

wiq) =f + a,Y.€ + ^2 ifr + 0{q'), (2.5) 

where ai and 0^2 are arbitrary constants. Here q"^ = qf + ql^ q = 2sin(g/2). 

Let us give some examples we will use in the following. The standard A^-vector model 
with hamiltonian 

H""'" = -NJ2(T^-(T^^^ (2.6) 
corresponds to w{q) = and thus we have ai = 02 = 0. Other possibilities are: 
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1. the Symanzik improved hamiltonian ||T5 



XfJ, 



4 



12 



' x+2fi I , 



for which we have 



1 



w 



(2.7) 



(2.8) 



and «! = 1/12 and a2 = 0; 
2. the "diagonal" hamiltonian |^| 



x+d 



(2.9) 



where d are the two diagonal vectors (1, ±1), for which we have 



w{q) = q - -qiq2. 





(2.10) 



and «! = 1/12 and a2 = —1/12; 
3. the perfect laplacian introduced in various works on the renormalization group R 



and recently revived in connection with the perfect actions [|I7 
1 



qfqi 



OO OO 

^ ^ 3^ ^ ;,?oo^.?oo + 27r/i)2 + (g2 + 27rhr {qi + 2nhnq2 + 2T,hy 

(2.11) 

for which ai = 1/12 and a2 = {n — 4) /(12/t). 

In general we will speak of tree-level improved hamiltoniansQ whenever ai = 1/12, 
a2 = 0: in this case, for 0, 



w{q) = q^ + 0{q''). 



(2.13) 



The hamiltonians (2/7) and (|2.11 ) for k = 4 are examples of tree- level improved hamil- 
tonians. 

In order to study the finite-size-scaling properties we must specify the geometry. We 
will consider here a square lattice of size L x T or a strip of width L with periodic 
boundary conditions in the finite direction(s). The large- limit of this model is well 

^Properly speaking we should speak of O(a^) tree-level improved hamiltonians. One can also consider 
0{a'^'') tree-level improved ones which are such that 



(2.12) 



for —> 0. We do not consider them here since tree-level improvement beyond O(a^) does not have any 
effect on the corrections to FSS at order For a perturbative study of this class of hamiltonians 

see Ref . [p5[ . Classically perfect hamiltonians are hamiltonians improved to all orders in a [p6[ . 
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known . The theory is parametrized by a mass parameter m| rp related to f3 by the 
gap equation 

where pi = 27Tni/L, p2 = 2Tm2/T and the sum extends over < < L — 1 and 
< ^2 < T — 1. The two-point isovector Green's function is then given by 

G.(.-.;L,r).K..„)^--E-^-^. (2.15) 

All other Green's functions are obtained from Gv{x\ L,T) using the factorization theo- 
rem 

In particular we will consider the isotensor (spin-two) two-point function 

GT{x-y;L,T)= ((o-, ■ o-,)^) - 1 = {a, ■ (Tyf + 0{1/N). (2.17) 



Beside the standard A^- vector model we will also discuss a mixed N-vectoi /RP 



N-l 



model ||2^, We will restrict our attention to nearest-neighbour interactions 



since only in this case the model is easily solvable in the large- iV limit. The hamiltonian 
is given by 



{l-r)cr^-(Ty + ^{(T^-cry) 



(2.18) 



where the sum is extended over all links (xy) and r is a free parameter varying between 
and 1. For r = we have the nearest-neighbour A^-vector model, while r = 1 corre- 
sponds to the RP^~^-mode\. Notice that for r = 1 the theory is invariant under local 
transformations (Tx — > €x(Tx, = ±1. Therefore for RP^~^ only isotensor observables 
are relevant. The limit we consider here corresponds to A^ — >■ cxd with r fixed. We 
mention that this is not the only case in which the model is solvable: a different large- A^ 
limit is considered in Ref. 



Also in this case the theory is parametrized by a mass parameter m\ rp related^ to j3 
by 11 

B= ^hAmlr? .219) 

where 

hAmlr) -^-^ , (2.20) 

with pi = 2niTT/L, p2 = 2n2Ti/T. 



is related to m^^T by Eq. ( |2.19| ) only for (3 > /3c (f"), where (3 = /3c (r) is a first-order transition 
line |4^. In the following we will be only interested in the limit /3 ^ oo, so that we will always use Eq. 



(Hi 
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The isovector Green's function is given by 

Gy{x - y; L, T) = — ^ -— . (2.21) 

^L,T[mL,T) ni,n2 P + ^L,T 

All other correlations are obtained using Eq. (|2.16|) . In the A^- vector case we can use 
the gap equation to substitute /l^tI'^lt) with [3. 

In this paper we will study the finite-size-scaling properties of various quantities. We 
define the vector and tensor susceptibilities 

Xv{L,T) = J2Gvix;L,T), (2.22) 

Xt{L,T) = Y.Gt{x-L,T). (2.23) 

X 

Using the explicit expressions for the two-point functions we get 



1 



XviL,T) = ^ ,^2^2 (2-24) 



Xt{L,T) - — — - — 2 — TT^TT^ / / \ , 2 — \2- (2.25) 

We want also to define a quantity behaving as a correlation length. In an infinite lattice 
there are essentially two possibilities. One can define the exponential correlation length 
from the large-x behaviour of a given two-point function^ G{x\oo): one considers a 
wall-wall correlation function 

G('")(l/;oo) = ^G((x,i/);oo), (2.26) 

X 

and then defines 

^(-p)(oo) = - lim , . (2.27) 

^ ' y^+oo logGW(?/; oo) ^ ' 

The mass gap yu(oo) is the inverse of ^'^^^^\oo). A second possibility is the second- 
moment correlation length ^^^\oo) that is defined by 

■^W(oo)]- (2.28) 
J 4 Y.xG[x]oo) 

The factor 1/4 has been introduced in order to have ^^^\oo) = ^^^^^-'(oo) for Gaussian 
models. 

We must now give the definitions in finite volume. Of course the exponential correla- 
tion can only be defined in a strip. In this case we can still use the definition (|2.27| ). For 
the second-moment correlation length we can use any definition that becomes equivalent 



•^Here and in the following wc will indicate the infinite-volume limit of an observable 0{L,T) with 
0(00). 
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to ( p.28| ) in the limit L,T ^ oo. We will consider here three different definitions: given 
a two-point function G{x; L, T) and its Fourier transform G{p; L, T) we define 



|^^(m,i)(2,,T) 



2pL 

-I 

87r2 ' 



1 



GiO;L,T) ^ 
G{po^; L, T) 
G{pox]L,T) 



+ 



my 

87r2 \ 



1 



G(0;L,T) 

. Li/2J Lr/2J 

--^ — E E 

4G(0;L,T) j=i_|^(i+i)/2j j=i-[(r+i)/2j 



G(0;L,T) y 
(z2+/)G((z,j);^,T) 



(2.29) 
(2.30) 



(2.31) 



where = (27r/L,0), = (0,27r/T) and [xj is the largest integer smaller than or 
equal to x. The third definition evidently coincides with ( p.28| ) for L,T —>■ oo. To verify 
the correctness of the other two definitions notice that Eq. ( p.28|) can be rewritten as 



Ml 



oo 



1 d"^ - 
TT^Gip; oo) 



(2.32) 



p=0 



Expanding in l/L^ it is easy to verify that both ^^"^■^-'(L, T) and ^'^'"'^^(L, T) converge 
to (|2]2|) for L,T ^ oo. Essentially (121291) and (|]30D are definitions in which one ap- 
proximates the second derivative of G with the difference at two nearby points. Thus 
these definitions converge to ^'^'"^(oo) as l/L^ (notice that G{x; L,T) — >• G{x, oo) expo- 
nentially). The third definition represents instead the "best" approximation on a finite 
lattice since converges to ^'•'"•'(oo) exponentially. This is indeed a general 

result that can be proved using the relation 



ihq)-fmmq), 



lL/2\ L-1 

E ^'m = 2 E 

j=l-[(L+l)/2j n=l 



valid for every function /. Here / is the Fourier transform of /, g = 27m/ L and 



(2.33) 



^ /27m\ 



cos ■ 



Tcn 

17 



for L even; 
for L odd. 



(2.34) 



If f{q) is meromorphic (as a function of the complex variable q) in the strip < Re g < 
271, periodic of period 27r, and with simple poles at q^, then we can evaluate this sum to 
obtain 

- /" (0) - T E . 2,-/,f. .r-,^M q,), (2.35) 
4 Y sm {qi/2) sm{Lqj2) 

where Ri is the residue of /(g) at q^. Thus the convergence rate is Lexp{—Lqj/2) where 
qj = miuj |Im(gj)|. For the specific case of the isovector correlation length one expects 
the nearest singularities (for /3 — * oo at least) to be at g = ±ifi{L) where /i(i^) is the 
mass gap. Thus we expect a convergence rate of Le'^^^^^^"^ . A general Green function 
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will not be in general a meromorphic function of q as cuts will appear as well. We expect 
however that the definition ( p.31|) will show the same exponential convergence rate. 
Using Eq. ( |2.33| ) we can rewrite Eq. ( |2.31| ) as 



Y: ^{G{q;L,T)-G{0;L,T 



2G{0; L,T) („^_„2)^(o,o) ^ 



X 



(2.36) 



Let us now give explicit expressions for the isovector correlation length: using the isovec- 
tor two-point function ( |2.15D we get on a finite lattice: 



1 ( W{P0^) ^ WjpOy) ^ 



POx 



Ply J 



m 



2 ' 
L,T 



+ 



Stt^ \w{pox) + 'm'^ 'w{pQy)+m?j 



(2.37) 
(2.38) 



rt2=l 



(2.39) 



In infinite volume we have instead: 



Am), \ 



(2.40) 



For the mass gap fiv{L) and the exponential correlation length C,v^^\l) we must solve 
the equation 



w{iiJ,v{L), 0) + m 



L,oo 



. 



(2.41) 



An explicit solution can be obtained only for the simplest w{p). For the hamiltonians 
we have considered in this section we have 




m 



L,oo 



1/2- 



1/2- 



for H'^'^ and H'^'^'s. 
for H^y"". 



(2.42) 

In our calculation we will only need the expression of fiv{L) for oo ^ 0. In this limit 
we obtain 



fJ,v{L) = rriL,, 



12 



(2.43) 



Isotensor observables are defined using the tensor two-point function ( p. 17] ). For the 
mass gap it is easy to verify that ^t{L) = 2jj,v{L). 
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3 A^-vector model 



3.1 The gap equation 

In this section we want to discuss the corrections to FSS for the hamiltonian 

Let us consider a fixed value of (5. Let and itll t be the mass parameters 
corresponding to (3 in infinite volume and in a box L x T. It is immediate to obtain a 
relation between rrioo and Indeed from the gap equation we obtain 

^lA^lt) = ^oo("^L)- (3-1) 

Now let us consider the finite-size-scaling limit, rrioo, ttll t — 0, L, T — ^ oo with mi TL = 
z and T/L = p fixed. Using the results ( [B.88| ) and ( |C.6| ) we obtain 

-fm{z;p)\l + logL + — I (3.2) 



m 



L,T 



with corrections of order 0(log L/L'^), where 
32 

Uz;p) =-e-^'^^«(^^''), (3.3) 

AmA^- p) = ^(12ai + 16a2 - 1) (32e-^'^^«(^^'') - z") , (3.4) 
Am,2(^; p) = 167r(12ai + 16a2 - l)Fo{z; p)e~^^Fo{z;p) 

-4(8ai + 8^2 - l)e-4"^o(^''') - Att (j^i(z; p) + 32e-'"^°(^''')Ai) . (3.5) 

The functions Fo{z] p) and Ti{z] p) are defined in the appendix, Eqs. ( |B.54| ) and (p.90| ). 
The function fm{z]p) is the FSS function for the ratio m^/m|^. As expected, it is 
universal (it does not depend on the explicit form of the coupling J{x)) and depends 
on the modular parameter p. The corrections instead are not universal. However the 
dependence on J{x) is very simple: the only relevant quantities are ai and 0:2 that are 
connected to the small-g^ behaviour of w{q) and Ai given by 



Ai 



+ 7^ aillpj + a2(p' 



2\2 



(3.6) 



The corrections behave in general as logL/L^, except when 12ai + 16a2 — 1 = 0. This 
cancellation happens for improved hamiltonians for which ai = and 0:2 = and 
also for many other hamiltonians that are not improved but nonetheless satisfy 12ai + 
16a2 — 1 = 0. To understand the relevance of this combination, let us introduce polar 
coordinates = qcosO, qy = qsinO. Then 

w{q) = + T^(12ai + 16^2 - l)?"" + ^(12ai - l)g'cos4^ + ©(g^). (3.7) 
lo 48 

Thus, if (12ai + 16^2 ~ 1) = 0, one cancels the first rotationally-invariant subleading 
operator, leaving a correction that is associated to a lattice operator that vanishes when 



10 



averaged over the angle 6. This last property is the reason why this quantity does 
not couple to the leading correction. This fact is not unexpected. Indeed the leading 
correction to scaling is usually associated to a rotationally-invariant operator (for a 
discussion for the two-point function in infinite volume see Ref. |p9|). 

For ai = and 0:2 = the expression for Am,2{z; p) simplifies drastically, becoming 



Am,2{z]p) = 47r(32e 



-47rFo(2;p) 



f 1 



V967r 



Ai 



Thus for improved hamiltonians there is the possibility of eliminating even the 
corrections choosing J{x) so that 

Notice that this condition is global, that is it does not only fix the small-g behaviour of 
w{q), but it depends on the behaviour of w{q) over all the Brillouin zone. 
A particular hamiltonian satisfying Eq. (|3.9|) is 



ffSym2 



Xfl 



+ 15a ) cr 



1 

12 



6a 



+ acr^ ■ <7 



(3.10) 



where a = 0.00836533968(1). The function w(g) is given by 



w(g) = + Ti^E^i + 



(3.11) 



For this hamiltonian the corrections to FSS behave as logL/L'^. Of course one could 



improve further: using a hamiltonian with w(q') = g + 0(g ) satisfying Eq. ( |3.9|) one 
could get rid also of the terms log L/L^; however the cancellation of the terms will 



again require a global condition of the type ( |3.9| ) . 

It is interesting to understand our results in terms of perturbation theory. Within 
the Symanzik improvement program the conditions ai = and 02 = are required 
for tree- level improvement: if the theory is tree- level improved , then the corrections 
behave as 1/L^ instead of logL/L^. In Ref. it was shown that the condition 

( |3.9| ) is necessary for improvement at one loop. The simplifying feature of the model is 
that, once the theory is one- loop improved, it is improved to all orders of perturbation 
theory. This explains why, if condition ( p.9|) is satisfied, corrections to scaling behave as 
logL/L^. As we shall discuss in the following section, for a generic model, for instance 
for a mixed 0{N)-RP^^^ theory, we expect only the term to be canceled so that 
the corrections to scaling would still behave as l/(L^logL). 

We have performed various checks of the expressions ( p.3| ), (|3.4|) and ( |3.5| ). First 
of all we have compared our results with previous work. For the strip fm{z', 00) was 
computed by Liischer |^1|. In this limp^oo ^i,i(^; p) = 0, using the exphcit 

expression for Fq{z;p), Eq. (p.54| ), and Eq. (|A.12| ) 



we get 



fm{z; 00) 



32 

— exp 



2n , TT^ 

2-fE + log — 

z 2 



exp 



-Aj^Koinz) 



(3.12) 



n=l 
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L 


Rexact{L] 2, 1) 


Rexpani^Lj 2, 1) 


Dexact{L; 2, 1) 


Dexpani^Lj 2, 1) 




4 


0.23892847 


0.24124682 


0.0649479 


0.0752812 


-0.428 


6 


0.23332246 


0.23379910 


0.0399609 


0.0420854 


-0.345 


8 


0.23028713 


0.23044079 


0.0264318 


0.0271167 


-0.303 


10 


0.22856980 


0.22863479 


0.0187774 


0.0190671 


-0.282 


12 


0.22751364 


0.22754600 


0.0140699 


0.0142141 


-0.270 


14 


0.22681775 


0.22683573 


0.0109682 


0.0110483 


-0.262 


16 


0.22633412 


0.22634493 


0.0088126 


0.0088607 


-0.255 


20 


0.22572118 


0.22572580 


0.0060806 


0.0061012 


-0.247 


32 


0.22497047 


0.22497124 


0.0027345 


0.0027380 


-0.232 


64 


0.22453987 


0.22453993 


0.0008153 


0.0008155 


-0.218 


128 


0.22441005 


0.22441006 


0.0002367 


0.0002367 


-0.208 


oo 


0.22435696 




-0.148 



Table 1: Values of Rexact{L; z, p), Rexpan{,L; z, p), Dexact{,L; z, p), Dexpan{,L; z, p) and 
(L; 2;, p) for the standard hamiltonian if'^*'^, p = 1 and z = 2. 



that agrees with the result of Ref. 

We have furthermore performed a detailed numerical check for the standard hamil- 
tonian i/**'^. Given L, z and p we have first computed m^rp = z'^/L'^, then /3 from the 
finite-volume gap equation P = XL^rijnLT) finally from [3 = Ioo{m1^): in this 
way we have obtained for each lattice size L and z the ratio Rexact{L; z, p) = ml^/m\ j,. 
Then we computed 

5i{L]z,p) = - — ^[Re^act{L;z,p) - Re^pan{L;z,p)], (3.13) 
log L 

where Rexpan{L; z,p) is the r.h.s of Eq. ( p.2|) . In this way we have tried to verify that 
indeed Si{L]z,p) at fixed z goes to a constant for L — > cxd. Numerically we find that 
6i{L; z,p) varies slowly with L and that the behaviour is compatible with the presence 
of 1 / log L and 1 / log^ L corrections. A better check can be obtained if we include the 
term of order log^ L/L^ that can be easily computed 

A„,3(^; P) = {12a, + 16a2 - 1)' [96e-s-^«(^^'') - 4z2e-4-i^o(.;p) + ^\ (3.14) 



Then we consider 
62(1; z,p) 



logL 



log^ L 

Rexact{L; Z, p) — Rexpan{L; Z, p) — fm{z] p)A„j^3(^; p) 



(3.15) 

In this case we should be able to verify that 

5,{L-z,p)^52,{z-p) + ^-^^^ (3.16) 

logL 
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for large values of L. The results for p = 1 and z = 2 (this value of z corresponds to the 
region where the corrections to FSS are larger) are shown in Table |l] where we also give 
the deviations from FSS, i.e. the quantity 



D{L- z,p) 



R{L]z,p) 

fm{z]p) 



(3.17) 



A plot of S2{L;2,1) versus 1/logL shows, as expected, a linear behaviour from which 
we can estimate 620(2; 1) ^ -0.148 and ^21 (2; 1) ^ -0.291. 

Let us now consider the limits z ^ 00 and z ^ 0. Asymptotic expansions of the 
FSS functions can be obtained using the expansions of Fq{z; p) and J-'i{z] p) reported in 
sections |B.2.1| and [B.2.2|. For large z we have 



fmiz;p) 

^m,l{z]p) 
^m,2iz;p) 



1 _ 2e-^/— - 2e-P\ — + Oiz'^/'e-',z-'/^e-n, 
^ z ^ pz 



1 



(12ai + 16^2 - l)z^ e 




27r 271 

z \ pz 



^0{^z^l\-\z^l\'f') 

T^'-- -- --2(^-z jZ , ^^pz I pz 



-(12ai + 12a2 - l)z' e-\ — + e-''\^] 
o \ V ZTT V zvr/ 

-2^ . , z^ ( , ^ /27r 

+ -(12a, + 16.2-1) log - e-y- + e-y- 



(3.18) 



(3.19) 



-0(^3/2g-.^^3/2g-p.^ 



(3.20) 



The result agrees with what is expected: for z ^ 00 the FSS function converges to 
1 exponentially ||61|. Also the corrections vanish in the same way and thus they are 
extremely tiny for large z. 

Let us now consider the perturbative limit (small z). For finite values of p, for 2; ^ 1 
and z <^ 1/p, we find 



fm{z\p) 

^m,l{z]p) 
^m,2{z;p) 



^r7(^p)V2'^-exp^ ^ ^'^'^ ■ -^^^^ 



pz- 



- — (12ai + 16^2 - 1) + O 

47ra2 



l-A7rz'Foi{p)+0{z^)), (3.21 



-47r/(p^2) 



-47r(l - 12«i)Fio(p) 



-Attz'' 



[1 - 12ai)Fn(p) + 2a2Foo(p) + ^ - Ai 

87r 



(3.22) 



Oiz^). (3.23) 



Here //(r) is Dedekind's //-function 162 



^TTiT/12 



n (1 - e'™") , 



n=l 



(3.24) 



13 



and Foi(p), -Fio(p) and -Fii(p) are defined in the appendix: see Eqs. (|B.62| ), ( |B.63|) and 

dm. 

For the strip the previous expansions are not valid. In this case, for small we get 

U^;oo) = ^e"2^-e-2-/Ml + ^C(3) + 0(^^)j, (3.25) 

A^^i{z; 00) = - j(12«i + 16a2 - 1) + 0(exp(-27r/2)), (3.26) 

2 

A^,2(^;oo) = — (12ai-l)-— (12ai + 12«2-l) + 0(22). (3.27) 

These formulae can also be used when p ^ 1. Indeed they approximate the FSS functions 
for 1/p < 2; < 1. 

It is interesting to obtain these expansions within perturbation theory (PT). The 
idea is to start from the perturbative prediction for rriooiP) 



m 



2 



32g4.Aog-4./3 ^ o(/3e-^"^)) , (3.28) 



and the perturbative expansion of mL,T{P)- 

z = L mL,T = 2^ • (3.29) 

n=l P 

Then we use the last equation to express perturbatively (3 in terms of z and finally we 
substitute the result in Eq. ( |3.28| ). In this way we obtain the expansions ( ^.21| ), ( pj.22D 



and (|3.23|) and the analogous ones on the strip. It must be noted that in this perturbative 
expansion the combination (12ai + 160:2 — 1) arises naturally: indeed it is the coefficient 
of the unique log L/L^ term that appears in the expansion. Thus 12ai + 16^2 — 1 = 
is the improvement condition of the renormalized perturbative expansion. 

To conclude the discussion we want to comment on the validity of PT: finite-volume 
PT is valid in the limit /? ^ 00 at L fixed while the FSS limit we are interested in 
corresponds to/3— >oo,L— s^ooat^; fixed. Thus our perturbative derivation of the 
FSS scaling functions involves an a priori unjustified extension of the validity of PT 
P^ , For the leading contribution this should be correct (naively because the result 
is L-independent) but the situation is less clear for the corrections: in this case the 
explicit calculation shows that the extension is valid also for the 1/L^ corrections, but, 
as we shall see in the next section, this is a special feature of the large- A^- vector model: 
in general the corrections to FSS computed in PT need a "resummation" to correctly 
describe the FSS regime. 

The functions Am,i{z] p) and Am,2{z; p) are reported in the figures |l], H, and | for 
the torus with p = 1 and for the strip p = 00 for the various hamiltonians we have 
introduced. 

From these plots one can immediately recognize a few basic facts: the region where 
the corrections to FSS are larger corresponds to 1 < 2 < 4 (the same has been found 
numerically in Monte Carlo simulations of H^^'^ with = 3 |]lOl)- In this region, for 
H''^'^ and H'^^"'^ and small values of L, say L = 10, the log L/L^ term gives a contribution 
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L 


Dexact{L] 2, 1) D 


3xpan(^Lj 2, 1) 


DexactiL; 2, 1) DexpaniL; 2, 1) 


4 


0.1363230 


0.1435180 


0.0039380 


0.0116907 


6 


0.0736526 


0.0752402 


0.0035162 


0.0051959 


8 


0.0463553 


0.0468952 


0.0023457 


0.0029227 


10 


0.0320544 


0.0322879 


0.0016198 


0.0018705 


12 


0.0235981 


0.0237157 


0.0011666 


0.0012990 


14 


0.0181629 


0.0182287 


0.0008802 


0.0009543 


16 


0.0144512 


0.0144909 


0.0006858 


0.0007307 


20 


0.0098297 


0.0098468 


0.0004483 


0.0004676 



Table 2: Deviations from FSS for H'^'^^ and H^y"": here p = 1 and z = 2. 



which is 2 — 4 times larger than the 1/L^ term and the corrections are positive. For these 
two hamiltonians the corrections become negative for large values of z (this can be easily 
checked from the asymptotic expansions ( p.l9| ) and ( p.20D ). They are also negative for 
j^std ^YiQ small-z region on the strip or on the torus for large values of p. Numerically 
we find that H'^^"'^ is the hamiltonian with the largest corrections, while H^^^ is the 
"best" one, as expected. 

We have finally checked that our expansion ( ^.21 ) describes well the corrections to 
FSS even for small values of L. In table |^ we give DexactiL; p) and DgxpaniL; z, p) for 
jjdiag ffSym p = 1 and z = 2. For the first hamiltonian there is good agreement 
even at L = 4, while for the latter there is a somewhat larger discrepancy, probably 
due to the larger spatial extent of the Symanzik hamiltonian. We have also computed 
DexactiL; p) for the same values of p and z for H^y^"^: for L = 4 (resp. 10) we get 
Dexact{L;'^A) = 0.0005743 (resp. 0.00001853). The corrections are estremely small (at 
L = 4 they are 100 times smaller than those present for if''*"'): improvement really 
works! 



3.2 Observables 

Let us now compute the FSS curves and the corresponding corrections for the observables 
we have introduced in Sec. |^. We will first consider the quantities that are obtained 
from the isovector correlation function, then we will discuss isotensor observables. 



3.2.1 Isovector sector 



From Eq. ( p.2|) it is immediate to obtain the finite-size-scaling curves and their leading 
corrections for the various observables. The susceptibility xv does not require any 
additional calculation since 

= ^ . (3.30) 

Xv{oo) mi^T 
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For the second- moment correlation lengths, neglecting terms of order log^ L/L^, we 
obtain: 



^(m,l) 




m 



L,T 



-^(ai + a2)— ^ 



(3.31) 



m 



l+3^(12ai + 12^2-1) 



(3.32) 



t(m,3) 



1 - 



+ 



4 ysinh2;/2 sinhp2;/2 



+ {Qp{L){z) + ^Qp(T){pz) 



(3.33) 



where p{M) is the parity of M (M = L, T) and 

47r2(l + p2) + 2pV 



^i(^;p) 
^2(^;p) 

Q even{z^ 
Qodd{z) 



{4n^ + z^){4tt^ + p^z^ 
+ 



(47r2 + ^2)2 
(12q;i + 12q;2 - 1 



(47r2 + p2;^2)2 ' 

z"^ cosh2;/2 



(12q;i + 12q;2 - 1)- 



192 sinh^ z/2 
cosh2;/2 



2 ^ - (4q;i + 4q;2 - 1) 



32 sinh z/2 



— («! + a2)- 



(3.34) 
(3.35) 
,(3.36) 
(3.37) 



192sinh^z/2 ' "^'8 sinh z/2' 

Let us now consider the asymptotic limit z ^ oo. In the FSS limit it is easy to obtain 



= 1 - 2e" 



Am,2) 



{L,T) 



OO 



j^(m,3) 



iL,T) 



FSS 
2 



FSS 
2 



FSS 




-2e 



-pz 




+ 0(z-^/^e- 



':Z' 



-3/2 -pz 



e-n. (3.38) 



= 1 



27r2(l + p2) 



(P'Z^ 



+ 0{z-% 



2 2 ^ 



-l/2„-p2 



)■ 



(3.39) 
(3.40) 



From these expansions one immediately sees that the FSS function for ^'"^^ goes to 1 
only as a power as z — > oo. The approach is very slow and indeed it reaches 1 at the 
1% level only for z ~ 60. This is extremely inconvenient for Monte Carlo applications: 
indeed in order to determine numerically the FSS curve one has to perform runs up to 
the value of z where the FSS curve becomes 1 within error bars: in this case runs with 
;2 pa 60 are required, which means that simulations on very large lattices are needed. The 
origin of these power corrections can be identified in the definition that approximates 
the infinite- volume ^y"^(oo) with corrections of order Xjl?: the Ijl? terms give rise to 
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the corrections of order Xjz^. The first definition should suffer from the same problem 
because also in this case ^'^^{L, T) converges to ^lp\oo) with corrections of order 
Instead Eq. (|3.38|) shows corrections of order 0{e~^ / y/z). This is a peculiarity due to 
the particular form of Gv{x) {Gv{x) is a free-field two-point function). However for 
different Green's functions terms of order l/z"^ are expected and indeed they are present 
for C,T^'^\ cf. Eq. ( p. 681) . As expected the FSS function for converges to 1 with 

corrections of order ze~^/^: in this case the FSS function is 1 at the 1% level already at 
z ^ 15. 

The large- 2; behaviour of the FSS-functions can be easily computed not only in the 
large- limit, but for all values of A^. The basic observation is that Gv{x; L, T) converges 
to Gvix; 00) with corrections of order LPe~^^^°°'^^ . Therefore, in order to compute the 
large-z expansion, one can simply replace T) with Gv{x;oo). The function 

Gv{x; 00) is well known in the critical limit. Indeed, if Gv{p', 00) is the corresponding 
Fourier transform, then, in the limit p —>■ 0, ^y^\oo) 00 with pC,l/^\oo) = Q fixed, we 



have H, H, g] 



G'y(0;oo) 



D{Q). 



Gvip; 00) 

The function D[Q) can be expanded in the limit Q ^ in powers of Q^: 



(3.41) 



D{Q) = E 



2n 



n=0 



(3.42) 



with hQ = hi = 1. This expansion converges up to the three-particle cut, i.e. for 
IQI < 3sm where Sm is defined by 



lim /iy (cx))^y"^ (00) 



(3.43) 



Sm is the ratio between the second-moment and the exponential correlation length. 
Moreover D{Q) has a zero in correspondence to the one-particle poles, Q = iis^- In a 
neighbourhood of these points, we have 



D{Q) = sj^ + l 



(3.44) 



Using these results it is straightforward to compute the FSS scaling curves in terms of 
y = L/^y^\oo) in the limit y 00. Disregarding terms of order y^e~^, ?/^e~^^ we obtain 



FSS 
2 



FSS 
2 



FSS 




2n 



2 L E:r=o&n(27r/y)2n +^y~^py\ 



1 - 



2s, 



(3.45) 
(3.46) 
(3.47) 
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In the large- iV limit 6„ = for n > 2, Sm = Sz = 1, and y = z m. the FSS limit, so that 
one recovers our previous results, Eqs. ( p.38|) , ( p.39|) , ( p.40|) . For generic values of 
numerical estimates of the various constants are reported in Ref. |]6^. The deviations 
from the large- iV values are extremely small: for = 3 one finds from a strong-coupling 
analysis H 62 = -1.2(2) ■ 10-^ Sm = 0.9994(1), = 1.0013(2), while a precise Monte 
Carlo simulation gives ||] = 0.9992(6). Using Eqs. (fctsp , (|3^ ), ( CT ), it is 



evident that the first definition is always the most convenient one except for extremely 
large values of y {y > 20 for = 3), where the deviations are extremely tiny. This 
is in agreement with the observation of Ref, 



7 ^ y ^ 10, had finite-size corrections larger than ^y'''''' . Using their data we can 

check the large-?/ behaviour of the FSS function of ^y^'^^ . We find that the data of Ref. 
they belong to the range 7 ^ y ^ 10 — are well described by the formula 



they found numerically that, for 

(m,l) 
V 

(m,3) 



35 



^(m,3) 



iL,T) 



l + a 



iL,T) 



exp 



2^{r'')(L,r), 



(3.48) 



a 



for N 
for 



(3.49) 



where 

1.023 ±0.012 
1.001 ±0.007 
in good agreement with our previous results. 

Let us now consider the corrections to scaling. The term proportional to log L/L^ 
is identical in all cases to Am,i{z]p), cf. Eq. (|3l^) . The contribution proportional to 
1/L^ depend instead on the definition of C,- In Figs. ^ and ^ we report the deviations 
from FSS for the three definitions for the standard and the Symanzik hamiltonians (A^ 1 
and A|_2 are defined in analogy with Eq. ( |3.2| )). Notice that for ^jj"'^-' and the 
corrections proportional to 1/L^ do not vanish even when ai = ^, a2 = and Eq. ( p.9|) 
are satisfied. This is expected since the second-moment correlation length is an off-shell 
quantity. Therefore the definition of the correlation length must be improved, as well 
as the hamiltonian. For instance, if one uses ^(T'^^ and the Symanzik hamiltonian one 
does not see any improvement: this definition has large corrections to scaling, and the 
behaviour is worse for the Symanzik hamiltonian than for the standard one. In this 
case there is a simple remedy to the problem: modify the definition in such a way that 



e(L,T) ^ U + OiL-\T-\L-^T- 



Analogously one could proceed for ^ 



(m,3) 



The 



second definition is automatically improved but this is a peculiarity of the large- A^ limit. 
Let us finally discuss the mass gap pv{L) and the exponential correlation length 
l/fiv{L). We have computed the FSS functions expressing them in terms of 
fiv{L) itself, i.e. using as variable x = fiv{L)L instead of z. We get 



^{exp) 



fn 



X] 00 J 



1± 



A^^i(3:; 00) 
L2 



logL 



A 



M.2 



x\ 00 



L2 



(3.50) 



where A„ Ax; 00) = A^ i(x; 00) and 



A^,2(a;;oo) = A^,2(a;;oo) 



^(12ai±12a2-l)^(x;oo) 



+-(12ai ± 12^2 - l)e-4"^«(^;~) . 
3 



(3.51) 
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Thus only the term differs from the expansion of im?^/m\rp. The asymptotic 

behaviour for large x is analogous to (p.20|) while for x — we have 



A„,.2(a;; oo) = — (12ai - 1) - —{I2a^ + 12^2 - 1) + 0{x^). (3.52) 
io o 



Notice that, since this quantity is defined on-shell, it is improved once the hamiltonian 
is improved. 

3.2.2 Isotensor sector 

Let us now consider the isotensor observables. The calculation of the FSS function for 
the isotensor susceptibility xt is straightforward. We obtain 



Xt(oo) V ^ 

with corrections of order log^L/L^ where 



L2 



^xT,i(^;p) 



647r dFo ^ 

z dz 



z\p)e 



-47rFo(2;;p) 



■(12ai + 16^2 - 1) 



In \ dz 



(3.53) 

(3.54) 
(3.55) 



A 



XT,2 



327r(12ai + 16a2 - l)i^o(2; p)e-^'^-^«(^'^) 
- 3)e-4-^«(^;^) - 47r. 

2567rAie-^"^°(^''') . 



4(28ai + 32^2 - 3)e-^"^«(^'^) - ^-nTxi^z- p) 
dJ^i/dz{z;p) 



(3.56) 



9Fo/9^(^;p) 

As before the logL/L^ corrections cancel if 12Q;i + 16a;2 — 1 = 0. The function A^^ 2(-2; p) 
simplifies considerably if ai = and 02 = . In this case 

\ -1" 



^XT,2(^;p) 



1 



967r 



Ai 



2567re"^"^«("''') - A-rz^ + 2^ 



9z 



>;p) 



(3.57) 



Therefore, if Eq. ( |3.9|) is satisfied, xt has only corrections of order log L/L^. It is 
straightforward to compute the expansions of the various FSS functions in the limit 
z ^ 00. Using the results of Sect. | B.2.1| and |B.2.2| we obtain 

1 + \/2t:z e 

-v2 



fxTiz;p) 
^XT,i(^; p) 



+ y^2^ e-P' + 0{z-^l'^e-\ z-^/^e-P'), (3.58) 



^XT,2(2;;p) 



— (12ai + 16a2 - 1) (^V2^ e"^ + ^27rpz e"^" 

+0{z^/^e-%z^/^e-n, 
__(12a, + 12a,-l)(^ 

v2 ^2 



(3.59) 



+ 



V^npz ) 



-j{12ai + 16a2 - 1) log ^ [e-'^2^z + e-P'^2^^ 

+0(^'/'e-^^'/'e-''"). (3.60) 
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As expected, f^^{z;p) behaves as z^e'^, but p differs from the value it assumes for 
other observables (see e.g. the large-z behaviour of xv)- Indeed, while the exponential 
behaviour is completely general, the power p depends on the observable. 

In Figs. H and |^ we report the graphs of A^^, 1(2; p) and A^^ 2(-2; p) for p = 1 and 
different hamiltonians. The behaviour is very similar to the behaviour of m| t/^'oo- The 
FSS corrections are quite small. The Symanzik hamiltonian shows the best behaviour, 
while the diagonal hamiltonian is the one with the largest deviations from FSS. 

Finally let us compute the FSS curve for the tensor correlation length. We will 
restrict the discussion to the standard action; the generalization to generic hamiltonians 
is straightforward but the final expressions are cumbersome. Moreover we will restrict 
our discussion to that is to the definition used in numerical simulations. We 

obtain 



ST 



00 



A 



L2 



logL + 



^gT,2(^;p) ' 

L2 



(3.61) 



where 



hri^'^p) = 192 Ei{z; p)e' 



^€t,i(^;p) 
^St,2(^;p) 

and 

Ei{z;p) = 

^2{z;p) = 



^3{z;p) 



1 



87r2 
1 



1287r3 



P 
1 



]_dFo 

2z dz 



'z]p) 



1 



+ 



1 



^3(^;p) 



+—^{z;p) 



Fsipz; i/p) 
1 



2z dz 
1 dFi 



Fi{z-p) 



+ 



F^{z-p) F,{pz;l/p)^ 



1 



87r2 
27r2 



2z dz 

OFq 

3 6z dz 
1 dFo 



z]p) 



p^Fi{pz-l/p)^ 
1 



I 1 



+ 



+ 



i^3(p^;i/p), 
1 



2z dz 



{^\P) 



^^'^\f^{z-p) ' p^F^{pz-Mp), 
' F,{z-p) F,{pz- l/p) J_ 
FKz-p)^ 



logp 



p2F|(p;.;l/p) levr p2F|(p^; l/p). 



(3.62) 
(3.63) 



A.,(.;p) + =4^ + 

A™,2(z; p) + - 4e-^'^^«(^^'') + 167rFo(^; p)e~^^^^^'-^p\ (3.64) 

-il^;; p) 



(3.65) 



(3.66) 



(3.67) 



Using Eqs. ( |B.57| ), ( [B.58D , ( |B.106| ), and ( |B.107[ ) it is straightforward to compute the 
large- 2; behaviour of the various FSS functions. We obtain 



f^A^'^p) = l-T^ — + ^::r^ — + 0[z 



15^2 



P 



63z^ 



P 



(3.6^ 



20 



^Ct,i(^;p) = -?^^ + T^^^(3-2p^ + 3p^)+0(z-^), (3.69) 




v2 



p^)^ + 0(log^/z4). (3.70) 



As expected f^^ [z] p) approaches one as 1/ z and thus it reaches the asymptotic value 
for z ^ oo within 1% only for z ^ 9. Notice moreover that /S.^^^2{z'-, p) diverges logarith- 
mically as z — ^ oo. This fact signals the non-uniformity of the expansion in z. This is 
not unexpected. Indeed, for each fixed z, we expect the expansion to be reliable only if 
'^L.T "C 1, i.e. if the correlation length is much larger than a lattice spacing. Therefore 
we expect the expansion to be valid only if 2; -C L. If 1 <^ z <^ L, log2;/L^ is a small 
number and thus the expansion is completely under control. 

The functions A^j,_i(2;;p) and Agj,^2(^;p) are reported for p = 1 in Figs. ^ and 
^ From these plots, comparing with the analogous graphs for other observables, one 
can immediately see that the corrections to FSS for ^j!^'^^ are quite large. This is 
particularly evident in the large- 2; region, where the log L/L^ term goes to a constant 
(for the isovector correlation length this term vanishes exponentially), while the 
term diverges as log z. 



4 Mixed 0{N)-RP^-^ model 

In this section we compute the FSS corrections for the hamiltonian ( ^.18[ ). As before we 
want to obtain a relation between and m\ at fixed /3. Using now the gap equation 
( |2.19| ) we have 



2 



(4.1) 



We will now discuss the FSS limit in which L,T ^ 00, (3 00, mi^T,'^oo ~^ with 
mL,TL = z and T/L = p fixed. At leading order we can disregard the terms ^ '"^"^ 
and m| 2"JL^r(^i r) denominators obtaining 



that implies iL^Tij^LT) = ^ooijnio)- Thus, at leading order, the relation between ttil^t 
and moo is identical to the one we have discussed in the previous section. Consequently 
the FSS functions for the hamiltonian (|2.18|) are identical to the FSS functions of models 



with hamiltonian (|2.1|) , as expected on the basis of universality^. The corrections will 

"^Under suitable assumptions (absence of vortices) that are verified in the large-A^ Hmit, one can 
prove that the FSS functions for the RP^^^ model with periodic boundary conditions are equal to the 
A'^- vector FSS functions with fluctuating periodic/antiperiodic boundary conditions In the large- 
limit the antiperiodic contribution vanishes, hence i?P°° has the same FSS functions of the iV-vector 
model. 
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however be different. Writing 



Ioo{ml^) = h,T{^L,T) + ^ "^i,T , (4-3) 

a simple computation gives 



^=o(z:t^-i) iTT^'rzjr^^- (4-4) 



2 \mi^T ) '^h,T{mL,T) - r 
Solving for fn'^/m'j^j, we get finally 



ml^ ^(^; P) [I + ^ogL + + 

where /m(-2; p), p) and Am,2(-2; p) are defined in Eqs. (|3.3|) , ( p.4|) and (|3.5|) , with 

«! = ^2 = and 

AL,i(^;p) = ^ (z^ - 32e-^'^^»(^^'')) , (4.6) 
7rr 



A:;,,2(^;p) = y(2Fo(;.;p)+r)(;.2_32,-4.Fo(.;p)^^ (47) 

A^{L;z;p) = (^^ _ 32^-4.^0 (.;p)\ L ^ . (4.8) 

™^ ' ''^^ 2 V y logL + 27rFo(;z;p) -vrr ^ ' 



The result ([4.5|) is quite different from Eq. (p.2|) . Indeed, while before the corrections 
had a very simple dependence on logL, now the corrections involve A^(L; z; p) that is 
not a simple polynomial in logL. Notice that, for large L at fixed z, ^"[^{L^ z\ p) behaves 
as 1/logL. Therefore, in the FSS limit, the corrections still behave as logL/L^. 

We should make a second remark about ^"[^{L] z\ p). It is easy to convince oneself 
from the asymptotic expansions, Eqs. ( |B.57| ), ( [B.59| ), and ( [B.69| ), that Fq{z; p) assumes 
any real value. Therefore, for each value of L, there is a value z^ such that the denomi- 
nator in Eq. ( [4. 8] ) vanishes, and therefore A^(L; 2; p) diverges. When L ^ 00, z^ ^ oo; 
more precisely, using Eq. (p.57| ), we have Zc ~ e~'^'''\/^L. This singularity is a signal of 



the fact that the expansion is not uniform in z. For each z the expansion is valid only 
when L ^ z, i.e. for m^^T <^ 1- In other words the expansion makes sense only when 
the correlation length is much larger than a lattice spacing. 

The corrections are larger for the mixed model than for the vector model. For instance 
Am p)/Am,i(-2; p) = 2r so that the logarithmic correction in the RP°° model (r = 1) 
is three times larger than the corresponding one in the A^- vector model (r = 0). For 
the values of L that are used in Monte Carlo simulations, say 8 < L < 128, however all 
terms in Eq. ( [4. 5] ) contribute to the FSS corrections. In table ^ we report for r = 1, 
p = 1 and z = 2, the same quantities reported in table |l[ Comparing the two tables 
we see that the corrections to FSS for RP°° are seven times larger than those for the 
A^- vector model in the same range of values of L. Only for L 50 (resp. L ^ 150) the 
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L 


Rexact^L] 2, 1) 


Rexpani^Lj 2, 1) 


Dexact{L] 2, 1) 


Dexpani^L^ 2, 1) 


8 


0.35709886 


0.34256990 


0.5916549 


0.5268967 


10 


0.29398737 


0.29002279 


0.3103555 


0.2926846 


12 


0.26809044 


0.26656021 


0.1949281 


0.1881076 


14 


0.25463070 


0.25391316 


0.1349356 


0.1317374 


16 


0.24665643 


0.24627448 


0.0993928 


0.0976904 


32 


0.22943124 


0.22941299 


0.0226170 


0.0225356 


64 


0.22561220 


0.22561112 


0.0055948 


0.0055900 


128 


0.22467784 


0.22467777 


0.0014302 


0.0014300 


oo 


0.22435696 





Table 3: Values of Rexact{L; z , p) , Rexpan{L] z, p), Dexact{L; z, p), and Dexpan{L; z, p) for 
the RP°° model (r = 1), p = 1 and 2 = 2. 



corrections are smaller than 1% (resp. 0.1%). To compare the corrections for the RP°° 
and the A^-vector model with the standard hamiltonian for all values of z, in Fig. ^ we 
report 



S{L; z;p) 



D. 



expan 



{L;z,p) 



r=l 



(4.9) 



Dexpan{L; Z, p)r=0 

for L = 128, 512 and p = I {D{L; z,p) is defined in Eq. ( p. 17] )). For these values of L, 
corrections are 5-10 times larger in the RP°° model. 

It is interesting to understand the origin of A^^lL; z; p) in terms of perturbation 
theory. First of all, for z small, (for simplicity we consider the strip case, analogous 
results are valid for generic p) we can expand, cf. Eq. ( p.69|) . 



^miL;z;p) 



2 TT + z log L 



1 - 



Tiz{2F. 



00 



TT + Z log L 



+ 0{z^) 



(4.10) 



This is not yet a perturbative expansion in z due to the presence of the term zlogL 
in the denominators. This term cannot be expanded in z since we are considering an 
asymptotic expansion at fixed z with L ^ 00. However, if we ignore this problem and 
expand the denominators in powers of z, we obtain 



^m(-^; ^1 P)pert 



Txr 



1 - - log L - z(2Foo - r) + 0{z^ log' L) 

TT 



(4.11) 



In general A5^(L; z] p) has a polynomial expansion in z with coefficients that are poly- 
nomials in logL: 

00 

(4.12) 



A;;,(L;^;p)p,,, = ^3^P„(logL)z", 

n=0 

where P„(x) is a polynomial of degree n. This expansion is clearly incorrect in the FSS 
limit L ^ 00 at fixed z. However Eq. ([4.12|) correctly describes the theory in the limit 
2; at fixed L. This is the limit in which PT works correctly [0, 4^, and indeed 
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the expansion ( 4.12 ) can be directly obtained with a perturbative calculation. Therefore 
our results show that in order to correctly compute the corrections to FSS one needs 
to resum the perturbative expansion. This reflects the fact that the perturbative limit 
z — > with L large and fixed does not commute with the FSS limit L oo with z fixed 
and small. It should also be noticed that the infinite series of logarithms appearing in 
Eq. ( [4.12 ) resums to give corrections of order l/(L^logL). This is a result far from 
obvious: in general one expects series of the form ( |4.12| ) to give powers of L, i.e. to 
resum to Lp^^\ where p{z) is some function of z. For = cxd no power is generated, 
but we have no proof that this will be true for generic values of A^. The only argument 
we have against the appearance of power corrections is based on a naive application of 
renormalization-group ideas. The corrections to FSS are due to the irrelevant operators 
of the theory. Since the A^- vector model is asymptotically free, operators have canonical 
scaling dimensions with logarithmic corrections. Therefore we always expect a behaviour 
of the form (logL)P/L^. 

Using the results of the previous section it is easy to obtain the FSS functions and 
their leading corrections for the various observables. For the isovector second-moment 



correlation lengths the expressions ( |3.31|) , ( p.32| ) and ( p.33| ) still hold with m^/m 



given by Eq. ( [4. 5] ) 

Xv{L,T) 



For the susceptibility xv we have 



Xv[oo) 



m 



1 + 



(32e 



-47tFo{z;p) _ ^2 



logL + 2nFo{z;p) 
log L + 2ttFo{z; p) - vrr 



(4.13) 



while for xt we have 



XTiL,T) 
Xt(oo) 



Finally for C,t^'^^ 



'xT{L,Ty 



\ Xt{oo 
r 

we have 



r=0 L 



1 + 



L2 



logL + 



(32e' 



-AirFo{z;p) 



L2 

logL + 27rFo(z;p) 
log L + 2-nFQ{z] p) - nr 



L2 



(4.14) 



ST 



iL,T) 



ST 



iL,T) 



(m,l), 



OO 



(m,l). 



OO 



1 + 



m,l 



L2 



logL + 



L2 



^miL;z;p) 
L2 



=0 



(4.15) 

To conclude our discussion we come back again to the RP°° model. In this case the 
FSS functions are usually reported in terms of s = L/^j^'^\l,T). Indeed, because of 
the gauge symmetry, one cannot define observables in the isovector sector. For any 
observable O, we define the FSS deviations 



Doms;p)=L' 



OiL,T) 



[ 0{oo) fa{s;p) 



1 



(4.16) 



where fo{s;p) is the FSS function of O expressed in terms of s. In Figs. |Ty and [11 
we report D^^^^L; z; 1) and D-^^{L] z\ 1) for L = 16 and L = 128. The corrections are 
extremely large if one compares them with the analogous results for the A^-vector model. 
This is especially true in the large-s region. Moreover the corrections are positive. These 



results are in qualitative agreement with the results of Ref. [47 
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5 Conclusions 



In this paper we have investigated the corrections to FSS in the large- hmit for a vast 
class of models and we have studied their relation with the improvement program of 
Symanzik. 

In the large- limit we find that the corrections behave as f {z,L) log L/L"^ where 
f{z,L) can be expanded in powers of 1/logL and is such that f{z, oo) is finite for all 
values of z. Thus, for large values of L, corrections behave as logL/L^. Tree-level 
improved hamiltonians have corrections to FSS behaving as the effect of the 

improvement is the cancellation of a logarithm. Subsequent perturbative improvement 
should give hamiltonians with corrections of order 1 / (L^ log' L) (I = 1 for one-loop 
improvement and so on). 

We have shown explicitly that the FSS limit and the perturbative limit commute in 
the calculations of the FSS functions but not for the calculation of the next-to-leading 
term. Corrections to FSS cannot be computed in perturbation theory unless an infinite 
series of logarithms is resummed. 

Finally we have investigated if there is any sign of large correction in the RP°° model. 
We find that this model shows deviations from FSS that are much larger than those of 



the A^- vector model. We find qualitative agreement with the results of Ref. |47 
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A Definitions 



A.l Functions G'fc(a), Hk{a) and 7ik{a] 
Let us define the following functions: 



(^2 ^ ^2^.-1/2 



E 

n=l 
oo -| 

V 



E 

?TI=0 



^2m^2fc-2m-l 



m 



5(1 



)fc+l/2 
1 



; (l-4n2)(n2 + a2)fe+i/2 



(A.l) 
(A.2) 
(A.3) 



The first and the third ones are defined for integers > while the second one is defined 
for k>l: for these values of k the sums converge for all values of a. The functions Gk{oi) 
and Hda) are known in statistical mechanics under the name of remnant functions (see 



App. D of Ref. m and Ref. 

We want now to discuss their behaviour for a — and a ^ oo. We will focus on 
those values of k that appear in our final results, i.e. to A; = 0, 1. 



25 



The expansion for a — is trivial. We obtain 



Gi{a 
Hi{a 
'Ho{a 

where we have used 



'2k 



g(-if^;jc(2.+i)^|) 



2k 



2E 

k=l 

oo 



(-1)^ /2A;^ 
k \ k , 



a 



^2k^ 



2k+2 



j:{-in2k+i)[^^ja2k+3) (^|) 

oo 

E(-i 



2k 



k=0 
oo 



'2k^ 
, k , 



l-21og2 + E^ 

° ^ — ' An 

n=l ^ 



4 5:(-i)^(2fc+i^ 

fc=0 



^2A;^ 
. k , 



a 



2k 



l-21og2 + V / ' 

^ ^ — ' An 



n=l 



n=l 



2q+l 



4^(1 -2 log 2)+ 4"^^^^''+^^ 

s=l 



(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 

(A.9) 



All series converge for |a| < 1. 

We want now to derive the asymptotic expansions for a oo. In order to obtain 
them let us derive a different representation for the functions Gfc(a) and Hk{a). 

Let us first consider Go{a). We rewrite it as 



Gr)(a) = lim 



''^' + '^ + 1(1/2 + 6/2)70 



(A.IO) 



Using the Poisson resummation formula [BH] we obtain 



Go{a) 



lim <^ -C(l + e) + 



a 



r(l/2 + e/2) V 2 



1 ^ + e 



r(l/2 + e/2) 



X 



2 

^exp 



n=l 



9 2 ^ 



X 



Taking the limit we get the representation 



a ]^ oo 

Go(a) = - log - - 7£ - — + 2 ^ Ko(27™a) 



(A.ll) 



(A.12) 



where Kq is a modified Bessel function |^^. The corresponding representations for Gi^ipi) 
and Hk{a), k > 1 can be obtained by integration and derivation of the previous relation. 
We obtain \M 



Gi(a) 



1_ 

12 + T 



a 1 1\ a ^ 1 



+ ^-lE- - ) - -J2 -Ki{2Tma] 



a/ ;;r^n 



HAa) = - 



1 An ^ , 

H nKi[2Tcna). 



oP- 20? a 



(A.13) 
(A. 14) 



n=l 
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This representation of the functions Gk{<y) and Hk{a) allows an immediate derivation 
of the asymptotic expansion for large values of a since 



KJx) 



1 + 



2x 



n 



4 



+ 0{x- 



Let us now consider the functions Tikia)- We rewrite 1-Lk{a) as 

1 1 



where 



-oo g— x'n^ 



E - 



Ait? 



(A.15) 



(A.16) 



(A.17) 



Now in the integral appearing in Eq. ( [A. 161 ) the relevant region for large corresponds 
to small values of x. Therefore we need the small-x expansion of g{x). First of all notice 
that 



0. 



(A. 18) 



n=l 



Then one immediately verifies that g{x) satisfies an equation of the form 

1 1 

g\x) + -g{x) - - E e— ^ = 0, 

n=—oo 

SO that, using Eq. ( |A.18| ), we get 

g{x) = -e-^l' / dy ey'' ^ e-^"^ 

^ -''^ n=-oo 

Using the Poisson resummation formula we can rewrite it as 

9{x) 



(A.19) 



(A.20) 



4 io ^ ^t^^ 

In the small-x region only the term with = is relevant so that 



(A.21) 



g[x) 



V^g-x/4 r ^ gS,/4 

Jo ^ 

' 2 S 4 

1 - 77 + t:^ - tttt: + 



6 60 840 15120 



+ 0{x^ 



(A.22) 



We obtain eventually for large values of a 



1 1 



+ 



TT 



1 



2«2fc+i 4 r(A;+ l/2)a2fc+2 



X 



k\ 



{k + l)\ ^ {k + 2)\ (A; + 3)! , (A; + 4)! 



6a^ 



60a4 



840^6 15120a8 



+ 0(a 



-2fe-12\ 



(A.23) 
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A. 2 The functions Mpq{z;p), Npq{p), Mpq{z; p) and Mpq{p) 
A second set of functions appear in our calculations. We define 

1 1 



Mpq{z;p) 

KM 



n=— oo 
oo 



4n2)(47r2n2 + z2)p/2 [exp(pV47r2n2 + z^) - 



y A i 

nP [exp(27rpn) - 1]?' 

oo I 

V - - 

(1 - 4n2)nP [exp(27rp?2) - l]? ' 



(A.24) 
(A.25) 
(A.26) 
(A.27) 



where g is a positive integer. We want to compute here the asymptotic expansion of 
Mpg{z; p) and A4pq{z; p) for z ^ and 2; — oo at fixed p. The first expansion is 
straightforward and we get for Mpg{z; p) 



Mpg{z; p) 



27r 



CO r) f) 



^ ^ (pz)" n,,.t^,=o ni\...ng\ 



+2z^ 



57r^ 477 



+ 0(^^). 



(A.28) 



For A4pq{z; p) the expansion is analogous with the substitution of Npg with Afpg. 
For large ;z let us consider first Mpg{z] p). We obtain 



Mpqiz;p)= 



(2vr) 



exp 



gpV47r2n2 + z^ + 0{e~^f""z^/^-P). (A.29) 



(47r2n2 + z2)p/2 

This last sum can be evaluated using the Poisson resummation formula [p6|. Define 



p 



(47r2t2 + ;,2)p/2 



exp 



-gpV47r2t2 + 



Then 



Mpq{z-p)^ J2 fC^Trn). 



(A.30) 



(A.31; 



For 2; ^ 00 we have 

f{uj) = ^{qpzy~P{uj' + 47r2g2p2)(2p-3)/4g^p 
Vz 

Therefore 



2tt 



Mpg{Z-p) 



1/2-p 



-qpz 



l + 0{z-') 



l + Oiz-^- 



(A.32) 
(A.33) 



/qp \2'n , 

For M.pg{z\ p) we will not need the explicit asymptotic behaviour. It is however easy to 
convince oneself that, for large z, Aipg{z; p) goes to zero faster than z~^e~''^^. 
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To conclude this subsection let us derive a set of relations for the functions Npq{p). 



First of all let us notice that Ni^i{p) can be related to Dedekind's ?7-function ||62 



riir) 



n (l - e^-^"") . 



n=l 



Indeed 



NiAp) 



oo oo -1 
n=l fe=l ^ 



log 



no 



-27rfep 



.A;=l 



7rp 



Following the same steps it is possible to prove that, for p > 0, 

(-l)P dP 



N-pAp) 



(27r)P dpP 



This relation, together with 



d 



-^Npgip) = -2TTq[Np^,Ap) + Np-i,,+iip)] 



{AM) 



(A.35) 



(A.36) 



(A.37) 



allows to prove that all Npq{p) with p < —1 can be expressed in terms of No^q>{p). From 
Eqs. ( |A.36| ) and (|A.37|) we obtain the following relations we will use in the following: 



n.,Ap) = NoAp) + NoAp)= ^ + 

N-2AP) + N-2AP) = N-iAp) + ^N^iAp) + 2iV-i,3(p) 



(A.38) 



1 d' 



\og7]{tp). 



Att^ dp 



(A.39) 



Finally, for p > 0, let us derive a relation between iV4p+3,i(p) and A^4p+3,i(l/p) that will 
allow us to compute explicitly A^4p+3,i(l). Let us start from 



E 



TT 2n 1 
- + 



a a e^'^° - 1 



(A.40) 



It follows 



E 



1 1 

+ 



-\a^P + 4) + ^C(4p + 3) + ^Ar4p+3,i(l/p) 

^C(4p + 4) + ^C(4p + 3) + -^N^p+^Ap) 



(A.41) 
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However, for integer p > 0, we can also compute the first sum as 

1 n^P+^ + {pmyP+^ 



E 



2p oo -| 

k=Q m,n=l yil^ll-P) 

E(-l)'^jiT^C(4p -2k + 2)C{2k + 2). (A.42) 
k=a P 



Comparing Eqs. (|A.41|) and (|A.42|) we obtain a relation between A^4p+3^i(p) and A^4p+3,i(l/p). 
For p = 1 we obtain 

1 1 1 2p 

iV4p+3,i(l)= — C(4p + 4)--C(4p+3) + — 5:(-l)\(4p-2fc + 2)C(2A; + 2). (A.43) 
27r 2 27r 

Particular cases are 

iV3,i(l) = ^ - ^C(3) ^ 0.001871373 , (A.44) 
197r'^ 1 

iV7,i(i) = Hsioo " 2^^^^ ■ ^^-^^^ 

Analogous relations can be gotten starting from the more general sum 

°° / 1 1 \ 1 

^ (pm)4p+2 j (n2 + p2m2)'?' ^"^'^^^ 

We leave the derivation to the reader. We will not need these relations here. 



B Asymptotic Expansions of Lattice Sums 

In this appendix we will study generic sums involving the lattice propagator for a Gaus- 
sian theory with arbitrary interactions. We will present an algorithmic procedure that 
allows to derive systematically the expansion in powers of 1/L of sums of this type. 
The results will be expressed in terms of the functions we have introduced in App. 0. 
The method presented in this Appendix applies to a square lattice L x T with periodic 
boundary conditions but generalizes easily to other types of boundary conditions. It 
can also be used to study lattice sums in more than two dimensions. In Sect. |B.1| we 
compute some preliminary one-dimensional sums; the general procedure is presented in 
Sect. Q 



B.l One-dimensional sums 
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B.1.1 The Euler-Mac Laurin formula 



The basic tool we will use is the Euler-Mac Laurin formula In its general form it is 
given by 



m— 1 
k=n 



dxf{x) - -{f{m) - f{n)) 



N r> 



.=1 m 
1 

'(2ivTT)! 



' dxp''+^\x)B,r,+i{x-[x\), 



(B.l) 



where are the Bernoulli numbers and Bk{x) are the Bernoulli polynomials defined 

by 13, g 

In the following we will be interested in sums of the form 

-| aL—l 



^ n=0 



(B.3) 



where p = 27m/L. We will try to compute the asymptotic expansion of the sum (|B.3|) 
for L —>■ oo with a fixed. It is easy to obtain from Eq. ([B.l|) the following formula: 



aL-l 



7 E fip) 

^ n=0 



/ //(p)-:t7(/(2^«)-/(0)) 
Jo 2,71 ZL 



{2N+1)\ 



with 









Jo 


Lp 


Lp ' 


2^ ~ 


2n 



lip), (B.4) 



(B.5) 



Thus, as long as the last integral in Eq. ( |B.4| ) is finite, i.e. /'•^^^^-'(p) is integrable in the 
interval [0,27ra], the previous formula gives the asymptotic expansion of the sum ([B.3D 
in powers of 1/L up to order An important case corresponds to a = 1 and f{p) 

periodic of period 27r, i.e. f{p + 27r) = f{p). In this case all the 1/L corrections vanish 
and we obtain 



E fip) = / //(p) + 0{L 



~2N-1\ 



(B.6) 



It is easy to prove that a similar result holds for generic (i-dimensional sums. If /(pi, . . . , pd) 
is a periodic function in all variables of period 27r and di^ . . . dd'^fip) is integrable in 
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[0, 2t{Y for all rzi, . . . , rzrf such that rii + . . . + Ud < 2N + 1, then, for L — >• oo with Li/L, 
L2/L . . . , Ld/ L fixed, we have 

E'e'--- E'/(Pi,...,P.)= / ^^^J{p)+0{L-'^~\ (B.7) 
where in the l.h.s. pi = 27mi/Li. 



B.1.2 Asymptotic expansions oij^n^ 

Here we will discuss the asymptotic expansion of sums of the form 



E 

n=l 



nP' 



(B.8) 



for L ^ 00. When p is a negative integer it is easy to perform the summation exactly. 
Indeed (g = — p) 



The simplest cases are 



E^^ 

n=l 



L 

E^ 

n=l 

E^' 

n=l 
L 



1 _ g"(i+i) 



1 - e° 



«=o 



E 

n=l 



-L(L + 1)(2L + 1), 



4 ^ ' 



(B.9) 

(B.IO) 

(B.ll) 
(B.12) 



Let us now consider the sum ( |B.8D with p > 1. Rewriting it as 

L 00 
^1 niZ^l 



[B.13) 



where C(n) is Riemann zeta function, and using the Euler-Mac Laurin formula for the 
second sum, we get the asymptotic expansion 



T.^,= av) - E §r(n + p - 1)L— . 



n=l 



(B.14) 



Finally, taking in the previous formula the limit p 1, we have the asymptotic expansion 



where 7^; ~ 0.577215665 is the Euler constant. 



(B.15) 
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B.1.3 Asymptotic expansions of + 

In this section we consider sums of tlie form 

L-l 



a 



2\fc+l/2 



[B.16) 



n=l 



Again we want to compute tlie asymptotic expansion for L — > oo witli a fixed. Let us 
first consider the case > — 1. In this case we rewrite the sum as 



L-l 

E 

n=l 



.2 ^ Ml/2 



y + 2™^2fc-2m+l 



k+1 

+ E 

m=0 L 



k + 1/2 \ 2m „2A:-2m+l 



m 



L-l 

la-- E 



n=l 



[B.17) 



We have already explained how to compute the last sums in the previous subsection. 
We will now discuss the first sum that we rewrite as 



Gk+i{a) - E 



n=L 



r^2 ^ Ml/2 



k+1 



, A; + 1/2V 2m^2fc-2m+i 



fB.18) 



where Gk{cy) is defined in Eq. ( [A.l| ). The last term appearing in Eq. (p.l8| ) can be 
easily computed using the Euler-Mac Laurin formula ( |B.1D . 

In the following we will need the previous sums for /c = — 1, 0. Explicitly we have 



L-l 



1 



Ev/^^^T^ = tL{L-l) + G,{a) + ^{\ogL + ^E)-^+0{L-'), (B.19) 



a 



L-l 



E 



1 



logL + + Go{a) - — 



For A; < — 1 the computation is straightforward as no subtraction is needed in this case. 
For /c = 1 we have 



L-l 

E 



(n2 + a2)3/2 
where Hi{a) is defined in Eq. ( |A.3| ). 



HA a) 



1 



1 



2L2 2L3 



(B.21) 



B.1.4 Computation of Y.{.P^ + a^)"'' 

In this section we will compute exactly sums of the form 

L-l 



E 



1 



for integer values of q. As usual, p = 2sin(p/2). 



fB.22) 
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If q is negative the summation is trivial as (here k>l) 



tf=[^iy. (B.23) 



n=0 

Let us now discuss the case g > 1. Consider first q — 1. Then 

L-l ^ ^ 



+ 2 5 ^ - cos(27rn/L) ' ^^"^^^ 



where uj — a /2. Then notice that 



hm / dz = 0, (B.25) 

r-^ooJdr cu — cos{27rz/ L) 

where Dn is the rectangle in the complex 2;-plane bounded by the lines z = —1/2, 
2; = L — 1/2, 2; = ±iR. Using the residue theorem we get 

Y — = , coth -arcchu; . (B.26) 

^0 a; - cos(27rn/L) x^J^^ V2 J ^ ' 

We thus end up with 

L-l 1 



^ -T^ n — — , coth 

^0 P + a q;a/4+7? 



L arcsh 



(B.27) 



Higher values of q can be handled by taking derivatives with respect to of the previous 
formula. 

B.1.5 Asymptotic expansion of Y^{p^ + m^)'^+^/^ 
Let us now consider sums of the form 

g(p2 + ^2)fe+l/2^ (B.28) 
n=0 

where p = 2 sin(p/2), p = 2iTn/L. We want to study these sums in the finite-sizc-scaling 
limit, i.e. for L — ^ 00, m — > 0, with mL = z fixed. To compute these asymptotic 
expansions we proceed in the following way. 

Assuming L even (the final result will not depend on this assumption) we rewrite 

J2{f + m'y^'^' = 2^f:'(p2 + m^)'^'^'+ (4 + m^)'=+i/^ + m^'=+^ (B.29) 

n=0 n=l 

Then let us consider the expansion of (p^ + z^ jl?)^^^^"^ powers of Ijl?: it can be 
written in the form 

34 



where ah{n^, z^) is a polynomial in v? and z^. Let us indicate with Rk,q{p, z; L) the sum 
of the first q terms in ( p.3CI|) . Then we rewrite 



m 



n=0 



k+1/2 



L/2-1 

2 E 

n=l 
L/2-1 



Z 



2 \ fc+1/2 



L2 



- Rk,q{p,z;L) 



+2 ^ Rk,gip,z;L) + (4 + 



2\fc+l/2 



2fe+l 



. (B.31) 



n=l 



We must then choose q. To fix its value we must decide the order in 1/L to which we 
want to compute the expansion. Then we fix q so that we can use the Euler-Mac Laurin 
formula for the first sum. It is trivial to reduce the second sum to a sum of terms of the 
form studied in the previous paragraph. 

We will now illustrate the method by computing the asymptotic expansion of 

L-l 



1 



V - ^ 



including terms of order l/L^. Since 
1 1 



p'^ + zyi'^ jp'^ + z'^/L^ 



1 + 



p 



24^2 + ^2/2,2 



+ 0{L- 



(B.32) 



[B.33) 



we rewrite 

L-l 



E 



Q Vp^ + m'^ 



L/2-1 

2 E 

n=0 
L/2-1 

2 E 



1 



1 



1 



p 



Vp^ + rn^ v^j?~+m^ 24 (p^ -|- m2)3/2 



+ — + 



1 



1 



L/2-1 

E 



p 



[B.34) 



Vp' + m2 12 (p2 + m2)3/2 74^^- 

The first sum can be computed up to order l/L^ using the Euler-Mac Laurin formula. 
We obtain 

L/2-1 r 



E 

n=0 L 



1 



1 



1 



p- 



yjp^ + Tii^ y/pP^^Trn^ 24 (jP' + m2)3/2 

2^11 



1 TT 7r 



4 96 167rL2 V 12 2 4 ^ 7r2 



TT 



Z 



1 



1 



1 



1 1 

1^2^ ' 48 2L2 V 16 ' 27r3 ' 167r 



TT / 1 1 
eZ W ~ 24 



B.35) 



The two remaining sums can be computed using the results of the previous subsection. 
We obtain finally 



1 



L 



L-l 

V - _ 

VP^ + 



TT 



log L + - log - + Go(-2/2vr) 



+ 



TT 



+^(^^i(^/2vr)-T7T 



z 



IOttL 



12/ 967r3L 

TT 4 



i^i(^/27r) 



log L + 7s - log - + -Go (2/27r) - - 
2 3 D 



+ 0{L- 



(B.36) 
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We will also need the expansion of ( B.29 ) for k = np to 0{L Using the same 
method we obtain 



L-l 



4L TT An 



n=0 TT 6L L L 



+ 2^ (logL + 7s - log - j + 0(L-3). (B.37) 
B.2 Two-dimensional sums 

In this Section we present our procedure to expand in powers of 1/L general sums with 
Gaussian propagators in the FSS limit. A general theorem for massless propagators was 



proved in Ref. |^. Here we will improve their result showing that only even powers 
of 1/L appear in the expansion and providing an algorithmic method to compute the 
various coefficients. 

B.2.1 Asymptotic expansions of J^Px^pf/'iP'^ + ^'^)~'^ 

In this section we present a general procedure to derive asymptotic expansions of sums 
of the form 

where Px = 27m^/L, py = 27my/T, the sum extends over < rix < L, < Uy < T, in 
the finite-size-scaling limit, i.e. for L,T 00, with T/L = p and mL = z fixed. 

First of all let us notice that rewriting p^^ = [{p'^ + m^) — p"^ — iv?Y "we can limit 
ourselves to consider sums with = 0, i.e. sums of the form 

1 p^^ 



The summation over Uy can be performed exactly using the results of section B.1.4 . It 
is easy to see that the result will be a sum of terms of the form 

L-l ^ r /I xi ^-c 

i2 _|_ /)oo2 



^ n,=0 ^ 



2Tarcsh \^\JPx + 'm' 



-1 , (B.40) 



for integers a, h and c > 0. If c is strictly positive it is simple to obtain an asymptotic 
expansion in powers of Indeed arcsh x = if and only if x = 0. Therefore, for 

L,T 00 the terms that contribute are those for which p ~ 0. Thus rewriting the 
previous sum as 

nx = lL/2\ \ / 

we expand the function g in powers of L at n, z, p fixed: 

9 [^^r^^pLj = Ej^^9m{n,z\p). (B.42) 
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The expansion of Eq. 



is simply given by 

oo 



- a—1 



E 

m=0 



1 



I2r. 



J2 9m{n,z^,p) 



(B.43) 



Let us now consider the case c = 0. If also 6 = we have discussed the asymptotic 
expansion in section [B.1.5| . Suppose now 6 7^ 0. Then define 



^ki = ^ / , 



ra=0 



n 



(B.44) 



and rewrite 



b/2 



nx=0 
L-1 



L 



E if. 



m 



a/2 



L-1 



+T E (pI + 



n:c=0 



m 



4 + + - Rbi 

a/2 

RhI- 



(B.45) 



Then choose / so that one can apply the Euler-Mac Laurin formula to the first sum: as 
the function is periodic of period 27i, as we observed at the end of section |B.1.1| (see 
formula ( [B.7] )) we can simply replace the sum with the corresponding integral. We thus 
obtain 



L-1 



m 



b/2 



1 



L-1 



tT.{pI + 



m 



a/2 



Rhi- 



(B.46) 



The computation of the remaining sums have been discussed in Section p. 1.5 . 
To illustrate the method let us consider a specific case, the sum 

/L,^("^^) = 77^ E -2] 2 - 

Using Eq. ( [B.27| ) we can perform the summation over Uy obtaining 

1 1 
E 

2 ^1 1 
+ T E 



(B.47) 



L nx=o J PI + m'^ J 4 + PI + m 



exp 



2Tarcsh ( '^yPx ~^ 



- 1 



(B.48) 
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The asymptotic expansion of the second sum is immediately computed: we get 



L-l 



Lj2 ; 

L n^=Q J PI + m'^ J 4 + PI + m'^ 



exp 



2T arcsh ( —yp^ + 'm'^ 



-1 



-) 



4 / ^ x2 



Ms,i{z-p)-2{—j Mi,i(z;p) -2M_i,i(2;;p) 



12L2 



2M_ 



1(2;; p) - 2 Mo,i(z;p) + f|^) V2,i(2;;p) 

4 



2nJ 
2 



27r/ 



+2M_2,2(^;p)-2f^) Mo,2(^;p)+ fT^) M2,2(^;p) 



2n 



+ 0{L' 



(B.49) 



Let us now consider the first sum. We want to compute its asymptotic expansion 
including terms of order l/L^. We rewrite it as 



Lj2^ ^ 

L n^=o Jpl + 711? JA+pl + m? 



■x 
L-l 

E 



1 L-l 



4 + p2+rn2 2 16 



L-l 



n,=0 a/K + ^ 



16L 



(B.50) 



The last two sums have been discussed in |B.1.5| . The first one, up to terms of 0(L~^), 
can be replaced by the corresponding integral. Expanding the integrand in powers of 
m? we get 



L-l 

E 



n,=o ^pl + n? 



4: + P1 + m? 
2 



1 1 

2 ^ 16 



/ -^2 2 



i-(l-log2)+ _^(l + 21og2)+ O(L-MogL). 



(B.51) 



It follows 



L-l 



L „^=o Jpl + m? JA + pl + m 



1 1 
2^ ^ 2^ 

TT 



log L + 7£; - log TT + - log 2 + Go ( — ) 



2ttJ 



6L2 



12 



2; 

2^ 



+ 



167rL2 



16L2 1927r3L2 
1 



log L + 7£; - log TT + - log 2 + -Go ( — 



0(L-^ log L). 

(B.52) 
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Using Eq. ( [B.49|) and the previous expression we obtain the following result: 



h,T{'m 
where 



2n 



\ogL+ Fo{z;p) 



,^^\ogL + ^F,{z-p)+ 0{L-HogL), (B.53) 



Fo(.; P) - i + ^ - ^ log ^ + Go {£) ) + ^M,,(.; p), (B.54) 



and 



TT / 1 



6 V12 
~16^ 



TT 



Z 



z 



z 



1 , TT 

- log h - 



16 1927r3 
2 2 



^1 — 
* 27r 





r 




~ 3 



\ 4 



M3,i(z;p)-2(;^) Mi,i(^;p)-2M_i,i(^;p) 



27r 
2 



2M_2,i(2 



;p)-2(^|-) Mo,i(2;;p)+ (^|-) M2,i(2:;p) 



\ 2 



+2M_2,2(^;p)-2(— j Mo,2(^;p)+ ( — ) M2,2(2;;p) 



27r 



. (B.55) 



Beside this sum we will also need 



TT L"^ 



^ - 2.G, ( ^ 
6 V27r 



— + {— 

4 167r3 ^ V27r 



47r 



1 2 



log L - log TT + - log 2 + 7^; - - + -Go — 



TT 
T2 



2M_ 



6 3 

\ 4 



2tx 



i,i(z;p) + 2f^) Mi,i(z;p)- f;^) M3,i(z;p) 



27ry 



27rV 
' L2 



2M_2,i(2;;p) - 2 f A")' Mo,i(z;p) + (^|-) V2,i(2;; p) 



2; 



27r/ 
2 



+2M_2,2(^;p)-2(— j Mo,2(2;;p)+ (^— j M2,2(z;p) 
12Fi(^;p)-|^) +0(L-4logL). 



1 3^2 1 



57r , 



+ 0(L-^ log L) 



(B.56) 



Finally we want to report the asymptotic expansions of Fq{z; p) and Fi{z; p) for 
z ^ and z +00. They are obtained using the asymptotic expansions of the 
functions Gk{z), Hi{z) and Mpq{z; p) reported in sections |A.1| and [A.2| . For z — >• +00, 
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we obtain 

-pz 





[z;p) = 


1 


log 


32 






32n 




'32 


Fi 


[z]p) = 





iTCZ 



(l + 0(.-)) + -^(l + 0(.-)), (B.57) 



y 24 72^ V ^ ^ 24 v/2^ V ^ ^ ^ 

(B.58) 

Let us now consider the perturbative limit. If p 7^ 00, we get for z <^ 1, z <^ 1/ p: 



1 



Fo{z;p) = + Foo{p) + z'Foi{p) + 0{z^), (B.59) 
Fi(^;p) = FM + z^Fn{p)+0{z'), (B.60) 



where 



i^oo(p) = ^(7i?-log7r + llog2)-^logr/(2p), (B.61) 

Foiip) = -T^C(3)-;^-^iV3,(p)-£,(iV,,(p) + iV,,(p)), (B.62) 

i^io(p) = ^- 3iV-i,i(p) + ^(iV-2,i(p)+iV-2.2(p)), (B.63) 

-^11 (p) = ^ ^7?; — logTT + i log2 — ^ + — logr7(ip) 

^'^^ 167r V 2 ^ 3/ 288 87r ^ 

2 

+ ^(A^o,i(p) + iVo,2(p)) - ^(iV_i,i(p) + 3iV_i,2(p) + 2iV_i,3(p)) 

= -1foo(p) - £Fio(p). (B.64) 

In the last formula we have used Eqs. ( |A.38| ) and (|A.39 ) that also show that -Fio(p) 
can be expressed in terms of derivatives of log ri{ip). For p = 1 we obtain the following 
numerical values: 



Foo(l) ^ 0.04876563317014130, 


(B.65) 


Foi(l) ^ -0.00386694659073721, 


(B.66) 


Fio(l) ^ -0.02924119479519021, 


(B.67) 


Fii(l) ^ -0.00376876379948390. 


(B.68) 


00) the previous expansions are not valid. 


In this case we write 



with 



i^o(^;oo) = ^ + Foo + z^Foi + Oiz^), (B.69) 
2z 

Fi{z;oo) = ^-^ + z'F,, + 0{z'), (B.70) 



Foo = ^(^7i?-log7r + ^log2), (B.71) 

Foi = -Y^C(3), (B.72) 

Fn = -Y^(7i^-log7r + ^log2-^). (B.73) 
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Finally let us comment on the duality property of the functions Fo{z; p) and Fi{z;p). 
The sum lL^j'{m'^) is clearly symmetric in L,T and thus it is a function ^{mL,T/ L, L) 
such that 

^mL, T/L, L) = ^{mT, L/T, T), (B.74) 

i.e. p, L) = ^{pz,l/ p, pL). This implies for the functions Fq{z] p) and Fi{z] p) the 
following relations: 



Mz;p) = —logp+Fo{pz;l/p), 
Zn 

z^ 1 
Fi{z;p) = ——logp+—Fi{pz;l/p). 

lOTT 



(B.75) 
(B.76) 



These equations provide a non trivial check for the correctness of our asymptotic ex- 
pansions and moreover imply the following relations on the expansion coefficients for 
z^Q 



Fooip) 
Foiip) 



1 



logp + Foo(l/p), 



Fio(p) = -Fio(l/p), 

p2 

Fiiip) = --^logp + Fn(l/p). 



(B.77) 
(B.78) 
(B.79) 

(B.80) 



The duality relation for Foo(p), -Fio(p) and -Fii(p) can be obtained directly from the 
inversion property of Dedekind's ?7-function 

r/(-l/r)2 = -zrr/(r)^ (B.81) 

where, in our case, we would identify r = ip. To prove directly Eq. (p.7^ ) one should 
use the relation obtained comparing Eq. ( |A.41| ) with Eq. ( |A.42| ) for p = 0. 



B.2.2 Asymptotic expansion of + 'Ti^) ^ 

We want to compute here the asymptotic expansion, including terms of order 0{L^^) 
of the sum 



1 



1 



— y - 



(B.82) 



in the FSS limit. Generic sums of the type ( |B.38D can be easily computed with the same 
technique. We assume that, for — vr < Pi < tt, w{p) vanishes only for p = and that in a 
neighbourhood of the origin w{p) has an expansion of the form (|2.5|) . Then we rewrite 



LT 



E 



1 



1 



w{p) + m'^ p^ + m? 
1 

LT 



+ 



tti EuPu + a2(F 



2\2 



E 



{p^ + im?Y 



LT ^ p^ + m? 



1 y^ C^lE^P^ + Q2(p^)^ 



[B.83) 
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Ao 


Ai 


ffSym 
}jSym2 


0.0322658881033520480 
-0.0471699346329274140 
-0.0564354728047190420 


0.00371978402668476 
0.00811339924292905 
0.00331572798724030 



Table 4: Values of Aq and Ai for various hamiltonians. 



Since we want to compute Xi r(m^) up to 0(L~^) we can substitute the first sum with 
the corresponding integral (cf. Eq. ( [B.7| )). Then, expanding the integrand in powers of 
m^, we obtain 

IlA^^) = Ao + tti 1 - - + 02 - T^Ai 



TT 



L2 



+lLAm ) - XT ^ if + m-)- ' ^^-^^^ 

where we have introduced 

= / (1^ - W + (^1? + cAff) ) . (B.86) 
and we have used the result (see Appendix C of Ref. |]38| ) 

^-^^ -1-1. (B.87) 



(27r)2 (p2)2 

We get eventually, using Eqs. ( |B.56| ) and ( |B.53| ) 
XL,T(m2)= --logL+ J-o(z;p)-^^(l -12ai- 16^2) log L + -^^1(2;; p), (B. 

where the neglected terms are of order 0(L~^ log L) and 

J^o{z;p) = Fo{z;p) + Ao, (B.89) 

J^i{z;p) = {l-l2ai)F,{z;p)+z^(^-Ai)+2a2Z^Foiz;p) + ^z'^{z;py, 

^^"^ ^ (B.90) 

Fo{z] p) and Fi{z; p) are defined in Eqs. (p.54| ) and (|B.55| ). Explicit values of Aq and 
Ai for the hamiltonians we have introduced in the text are reported in table ^ 

Using the expansions of Fq{z\ p) and Fi{z; p) (see previous section) we can easily 
obtain the asymptotic expansions of J^i{z; p). For large z we obtain 

2 22 

^i(z;p) = -(12ai + 16a2-l)7^1og|---^(8ai + 8a2-l)--2'Ai 



.l(12a, + 12a. - D- + ^) + 0(.'/V^.3/^e-). 



B.91) 
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For finite p and z <^ 1, 2; <^ l/p, neglecting terms of order z*^, we liave 



p) = (1 - 12ai)Fio(p) + y + 
wliile on tlie strip, for 2; — > 0, we obtain 



;i-12ai)Fn(p) + 2a2i"oo(p) 



^•1(2;; 00) = (1 - 12ai) + — (12ai + 12«2 - 1) 
72 io 



:i - 12ai)Fn + 2^2^00 + ^ " Ai 



57r 



Ai 



(B.92) 



0{z^). (B.93) 



B.2.3 Sum for the tensor correlation length 

In this section we describe the computation, in the FSS limit, of 



E 



1 



-^^'^v [{px - Po) + Pi + m'^][pl + Pi + m'^ 



where po = 27r/L and, as before, p^ = 2Tmx/ L and Py = 2TTny/T. 
First of all we rewrite Eq. (|B.94|) as 



then we sum over py to get 



9 L-l 



E 



1 



1 



^ ^0 {p^pof - p^ VP'^ + mW4: + P^ + m'^ 



exp[2Tarcsh(|v^p^ + m^)] — 1 



(B.94) 



(B.95) 



(B.96) 



where we have simplified the notation using p instead of px- 

The contribution due to the second term in curly brackets is obtained by simply 
expanding in powers of (see the discussion of Eq. ( p.4CI| ) for c 7^ 0). The remaining 
term requires more care. Assuming L even, we rewrite 

E 



L ^0 Tp^f _ p2 ^p' + m'^A+p' + 



LpQ a/4 + m2 \m 1/8 + 



L/2-1 



1 - pV2 



L (po — 4 sin^ p) ^p'^ + m^-\/4 + p^ + m' 



7 ■ 



[B.97) 



Consider now the last sum and notice that, for L ^ 00, m ^ 0, beside the singularity 
at p = 0, there is an additional singularity at p = vr. Using the fact that 



L/2-1 



E 



\ pQ — A siv? p 
43 



1 



fB.98) 



and keeping only those contributions that do not vanish for L — >■ oo we rewrite Eq. 
(IRQTI) as 



1 1 1 

+ 



L/2-1 



+ 



1 - pV2 



L -4sin2p V4V2 ' ^p'^ + m^^A+p 



[B.99) 



In this way we have removed the singularity for p = vr. The remaining part of the 
calculation follows the lines we have presented for Eq. (|B.40D when c = 0. We subtract 
to the sum the first two terms of the asymptotic expansion in and then replace the 
sum with the integral. Explicitly, if we define 



we obtain 



2(p2_ V) 



^ + ^--. ^ + 



12 4 6 24p2 + rn2 



+ 



2^2 



00) 



1 



+ 



1 1 



zpl V4 + m? 2V2Lpl 
1 

/o 27r 



/■'^ dp 

+4 / 7e- 



1 ^ 1-pV2 



4sin^p \4v/2 pa/4 + 
+ 7 E ^(p,Po,m2)+ 0(1/L). 

n=l 



i?(p,0,0) 
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The last sum can be dealt with following the strategy of section |B.1.3| . We get finally 
2 ^1 1 1 ^ 
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TT 
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r 2 oo 
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1 2 

log L + 7£; - log TT + - log 2 - - 
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1 oo r 1 



27r2 ^1 i (1 - 4r22)V47r%2 + ^2 ^ 3 



2^2 



1 



vr n 



6 24 47r2n2 + 



-]■ 

8nj 
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Collecting everything together and introducing the functions TCk(a) and Aipg{z; p) de- 
fined in sections |A.1| and |A.2| we have 



/2,LT(m2) = L^Fsiz; p) + — logL + F4(;2; p) + OOogL/L^), 



where 



IGtt 



1 1 

+ 



47r22; 4:71^ 



Ho 

44 



27r 
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(B.103) 



(B.104) 



and 
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1 1 
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7B-log7r + -log2-- 
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247r V27ry 
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247r 
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277 



+2A1_2,2(2:;p) -2 f;^yA<o,2(^;p) + { M2,2{z\ p) 



2TiJ 



2-n 



. (B.105) 



To conclude this section we give the asymptotic expansions of ^3(2;; p) and F/^{z\ p) for 
large and small values of z. The necessary formulae for the derivations are reported in 
sections |A.1| and |A.2|. For large z we have 
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27r2 
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327r 
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0(z-«). (B.107) 



where 



Fsoip) 
FM 



27r2pz2 

77^ + i^4o(p) + 0(2;2), 
bpz'^ 



\(l-21og2) + -^ + -i,M,i(p), 



47r3' ' 247r2 2ti^ 

^ -(7^-log7r)--^(2 + log2)+ ^ 
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167r 



967r 



72 247r2p 



+^{N,M+^iM - 2Ar_i,i(p)) + ^(Ar_2,i(p) + Ar_2,2(p)). (B.Ul) 



On the strip, for small z, we have 
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Fa{z\ 00) 
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+ 0{z^). (B.113) 



45 



C Asymptotic expansion of lattice integrals 

In this section we want to discuss the asymptotic expansion for mg ^ of the integrals 

d^p 1 



(27r)2p2 _|_ ^2' 
d^p 1 



(C.l) 
(C.2) 



(27r)2 w{j)) + mo' 

More general integrals can be discussed following the same method and using the results 



of App. C of Ref. The expansion of /oo(^o) easily obtained from its expression 
in terms of elliptic integrals p7| 



2 ' -.k' ' 



TT 4 + mQ y4 + mQ j 



(C.3) 



To obtain the expansion of X^im^ let us proceed as in section |B.2.2| . We rewrite 



2oo("^o) 



dp 
(27r)2 

+ 



1 ^aiE^P^ + a2(p')' 



W(p)+m2 p2 _|_ ^2 (^p2 _|_ j.^2^5 



dp 



(27r)2]32 + m2 7 (27r)2 (p^ -|- 77^2^^ 



2\2 



(C.4) 



If we want to compute the expansion neglecting terms of order O(mglogmQ) we can 
expand the first integral in powers of mg. Then, using (see App. C of Ref. |38]) 



dp E^P^ 
(27r)2 (p2 + rnlf 



l-- + ^(31og^ + 2) + 0(m^logmg) 



(C.5) 



we obtain 



+^(1 - 12ai - 16a,) log ^ + ^(1 - 8ai - Sa^) - m^A^, (C.6) 



2 2 



where Aq and Ai are defined in Eqs. (|B.85| ) and (|B.86| ) 
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Figure 1: Am,i(-2;p) for the standard hamiltonian H^*'^ for p = 1 and p = oo. The 
dashed hnes correspond to the asymptotic expansions ( p.l9| ) and ( p.22D . A.^ i {z; p) for 
}jdiag jg obtained by multiplying the vertical scale by 4/3. 
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Figure 4: A^^^2{z; 1) for H^^'^ for the three different definitions of second-moment corre- 
lation length. 
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Figure 5: A^^, 2(-2; 1) for H^^"^ for the three different definitions of second-moment cor- 
relation length. 
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Figure 7: A^^ 2(-2; 1) ioT H^^'^, H^"^"^ and H'^^"-^. The dashed hnes are the large-z asymp- 
totic expansions, Eq. ( p.60| ). 
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Figure 9: S{L; z; 1) for L — 128 and L — 512. For L ^ oo, S{L; z; 1) converges to 3 for 
all zj^Q. 
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